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1. Introduction 

The calculation of order parameters in solvable statistical models has by now a 
rather long history. It began with Onsager's result for the bulk magnetization of the 
square lattice Ising model [1, 2]. Quite some time later, McCoy and Wu obtained the 
order at a (10) surface of the same system [3, 4], by which the general study of surface 
critical behaviour began. Again it took some time until it was realised that this 
behaviour depends on the shape of the surface and that an interesting phenomenon, 
namely continuously varying exponents, appears at corners [5, 6]. Thus Ising models 
with corners have been studied over the last 12 years [7] but so far there has been only 
one completely analytical calculation, namely for a square lattice with a 90° corner 
formed by two (10) surfaces [8]. 

The calculations in all these cases can be (and have been) done using row transfer 
matrices. One then has to find, for example, the matrix element of a spin operator 
between two asymptotically degenerate eigenvectors of this matrix. However, while 
this is rather simple for the (10) surface [9, 10], especially in the anisotropic (Hamil- 
tonian) limit [11], it involves the solution of an integral equation for the bulk and the 
corner problem [10]. Thus an alternative and possibly simpler method is desirable. 

Such a technique has become available in the form of Baxter's corner transfer 
matrix (CTM) [12] . In fact, the bulk order parameters in most of the more complicated 
solvable models have been obtained through it. In this method, a square lattice is 
divided into quadrants and the CTM is the partition function of one such quadrant 
with the variables at the outer boundary held fixed. One then can calculate directly 
the expectation value of the central spin in the form of a trace. A basic feature is the 
simple structure of the CTM eigenvalue spectrum in the thermodynamic limit [12] 
which leads to simple infinite product formulae for the order parameters [13]. 

While in the early CTM calculations the spectrum was used directly a somewhat 
different approach has been developed recently by the Kyoto group [14-16]. In 
this approach, so-called vertex operators play a fundamental role. These are half- 
infinite lines of vertices in vertex models, or lines of couplings in an Ising model, and 
hence can be viewed as half-infinite row transfer matrices. If they differ from the 
lattice couplings only in their anisotropy (spectral parameter) , then they have simple 
commutation relations with each other and with the CTMs. This can be used to 
"rotate" such lines in a lattice and to obtain in this way functional equations (also 
known as g-difference equations) for order parameters or even correlation functions. 
(Actually the idea of moving lines about in a general way was introduced by Baxter 
for the eight vertex model [17], a concept called "Z-invariance".) Such equations can 
be solved in some cases in a proper elliptic functions parametrization and this leads 
to the same type of infinite product formulae as are observed in cases where direct 
calculation is possible. 

So far, this technique has been used mainly for obtaining some bulk order pa- 
rameters, e.g. in ABF models [18,19]. Only recently, vertex models with straight 
surfaces have been considered, and the spontaneous magnetization at a free end of 
the XYZ spin chain has been obtained [20] . An additional ingredient needed here are 
boundary Yang-Baxter equations and the corresponding K-matrices. 
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In this paper we want to consider the square lattice Ising model and show how, 
using this approach, the various types of order parameters can all be obtained in a 
unified way. We recover all previously known results and derive several new ones. 
Although the Ising model is a special case of the eight- vertex or the ABF models, we 
think that it is useful and warranted to treat it completely on its own. 

Thus we will first outline the technique of forming the functional equations in a 
spin formulation and show how to recover once more the Onsager formula. Then we 
describe how the procedure has to be modified for a system with a free (11) surface, 

1. e., one cut diagonally with respect to the couplings. It turns out that, instead of 
complicated Yang-Baxter boundary equations, one needs only a simple elementary 
relation here. The structure of the functional equations remains essentially the same 
as in the bulk and we obtain from it the magnetization in the first and second row. 
We then apply the procedure to a 90° corner formed by two (11) surfaces and obtain 
the magnetization at three different corner points. The solution for this corner is 
especially interesting in the critical region, because the corner exponent (3c depends 
on the anisotropy. We obtain this dependence, which previously was inferred by 
rescaling onto an isotropic system [5 - 7] , directly from our analytical formulae. 

Finally we show that the method also works for a lattice with a (11) surface and 
a layered structure of the couplings perpendicular to it. An appropriate choice of the 
layering then produces (10) surfaces, so that these can also be treated. In particular, 
one can obtain the other type of 90° corner in this way and derive the corresponding 
magnetization formula. This completes our unified treatment. A concluding section 
contains some additional remarks and an outlook. 

2. Basic constructions 

2.1 Parametrisation. We consider an Ising model defined on a regular square lattice 
taken on the diagonal, and parametrise the couplings using multiplicative rapidity 
variables d, (^2 attached to lines running in the horizontal and vertical directions. 
Except in section 6, we do not account for the horizontal and vertical rapidities 
individually, simply writing ( — Ci/C2- 



There are two types of coupling as shown in figure 2.1. Our parametrisation uses 
the standard formulae [12] 




Figure 2.1: Ising couplings K and L. 



sinh2iir = — zsn(m), 
cosh2i^ = cn(m). 



sinh 2L — —i sn{il' — iu) — i ns(m) / k, 
cosh2L = cn(i/' — iu) = ids{iu)/k. 



(2.1) 
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where sn, cn are the usual Jacobian eUiptic functions, of modulus 

k = l/(sinh2ii:sinh2L). 

This parameter, already introduced by Onsager, measures the temperature. The 
quarter periods of the elliptic functions are / and /', where / = I{k) is the elliptic 
integral of the first kind and /' = I{k') with k' = y/l — k'^. The anisotropy of the 
couplings is measured by the variable tt, or equivalently by 

C, = exp(— 7ru/2/). 

Furthermore we will use the quantities 

q = eM-^l'/I), ^ = 

In these variables the crossing symmetry (rotation through 90°) L{u) = K{I' — u) 
becomes L{Q = K{x/C). In our calculations we shall particularly need an expression 
for tanhL(C) in terms of infinite products. It is 

i-q'/yC; q'){-q'/H; q'WVC; q'){q'^K; q') ' ^ ^ ^ 

where {z;p) is the infinite product 

oo 

{z;p)^l[{l-zp-). 

n=0 

The identification (2.2) follows from standard elliptic function identities, together 
with the fact that the Jacobi theta function H(u) has the infinite product expansion 
H{iu) = iq^^^z~^/'^{q'^; q'^){z; q'^){q'^/z; g^), where z = exp(— ttw//). We note for later 
use the relation 

""^^^^ -fca/2g^^, (^;g^)(gV;;g\ ^ . = exp(-W/). (2.3) 



H{iu + I) cn(iu) {—z',q'^){—q^/z;q^) 

2.2 Star-triangle relation. Integr ability depends on special properties of the Boltz- 
mann weights, particularly the star-triangle relation shown in figure 2.2. The filled 
circle represents a spin which is summed out, while the scalar factor is determined by 
the normalisation of the Boltzmann weights. The condition on the rapidities in the 
three couplings is given by Baxter [12]; taking into account the present notation it is 





Figure 2.2: Star-triangle relation. 
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Cs = C1/C2. (2.4) 

This relation is important in choosing appropriate couphngs when deriving func- 
tional equations. For vertex operators, we shall use the star-triangle relation in an 
iterated form shown in figure 2.3. From figure 2.2 we see that H is of type L, M of 
type K; this is discussed further in section 2.4. Notice that in each pair of operations 
the individual scalar factors R are mutually reciprocal, and therefore cancel. 




K2 L2 K2 L2 K2 K2 L2 K2 Lj Kj K2 Lj Kj Lj Kj 



Figure 2.3: Iterated star-triangle relations. 

2.3 Corner transfer matrices. A CTM is defined to be the partition function of a 
quadrant of the lattice, as shown in figure 2.4. The direction of transfer (action as 
an operator) is shown by an arrow. There are two sectors, Neveu-Schwarz {NS) and 
Ramond (i?) according as the centre of the lattice is taken at a single spin site or not; 
thus A^^^\Q is diagonal with respect to the central spin e. We fix the boundary 
spins to +. 




Figure 2.4: Corner transfer matrices for the A^^" and R sectors. 

The really remarkable properties of the CTMs emerge in the thermodynamic 
limit. In this limit, we expect that the entries of the CTMs will scale as the partition 
function per coupling raised to the power of the number of couplings, and we assume 
that this factor is divided out to give a convenient finite limit. 
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These normahsed infinite CTMs are generated by spin chain operators H^^^^ 
and H^^')-. 

Prom this it is evident that the CTMs are symmetric and that they satisfy the group 
property 

^(^^)(Ci)A(^^)(C2) = A(^^)(CiC2), A(^\Ci)A^''\C2) = A(«)(CiC2). (2.5) 

The eigenvectors of H^^^^ and H^^'^ span the Neveu-Schwarz and Ramond sector, 
respectively. Moreover, the Neveu-Schwarz sector divides into two subsectors labelled 
by the value of the central spin e. Since it is straightforward to diagonalise H^^^^ 
and H(^^ using fermions [21, 22], analytic formulae may be given using properties of 
elliptic functions. In particular, the eigenvalue spectra are equally spaced, although 
we do not make explicit use of these facts in this paper. 

Finally we note that one may define CTMs for the other three quadrants of the 
lattice, North-East, South-East, South- West. They are simply related to the one we 
have defined (North- West) by 

^NE (C) = ^NW i^/O > ^SE (C) = ^NW (0 , ^SW (C) = ^NW i^/O ■ 

provided that ( retains its meaning as the ratio of horizontal and vertical rapidities. 
The prime denotes the transpose (which reverses the direction of transfer); this is 
the form in which we construct ^ne from ^nw in the non-symmetric case needed in 
section 6. 

2.4 Vertex operators. The central technique of this paper is to employ "vertex oper- 
ators" which map between (intertwine) the Neveu-Schwarz and Ramond sectors. We 
are concerned with the physical interpretation of vertex operators as semi-infinite row 
transfer matrices V^{Q and W^{Q, and in the functional relations which they satisfy. 
The two types of operators are shown in figure 2.5. Since we shall often rotate these 
operators to other orientations, it is important to note that they retain their under- 
lying definition as Ve{K, L), We{K, L), from which they derive their dependence on ^ 
via (2.1). One sees immediately from the crossing symmetry that the vertex operators 
satisfy 

K'(C) = w-e(x/C), w^iO^v^ix/c). 




Figure 2.5: Vertex operators V'e(C) ^nd W^{Q. 



Consider now the left hand side of figure 2.6, which shows the composition of 
the action of A^^^\Q with Ve(^). For a finite system this simply gives a Ramond 
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corner, but in the infinite case we must be more careful; A^^^\() maps the Neveu- 
Schwarz sector into itself, after which V^(C) maps it to the Ramond sector. So, on 
the right-hand side of the equation, we have inserted a trivial operator which makes 
the same mapping just by re-labelling the spins (the dashed lines indicate simple 
identification). However this operator is the anisotropic limit of V^{Q at C = 1, so we 
obtain the functional relation 

V,iOA(''^\0=A(''\C)Ve{l). 

By eliminating Ve(l) using the group property of CTMs, we obtain the general "boost 
property" 

^e(CiC2)^(^^nCi) = A^''HCi)VeiC2). (2.6a) 

This boost property of CTMs has been the subject of some investigation [23, 24], but 
its real utility only emerged recently [14-16]. 




Figure 2.6: Mapping from Neveu-Schwarz to Ramond sectors using K(C)- 

One can write a similar relation for the second type of vertex operator, which we 
shall denote We{C). The picture is shown in figure 2.7, and an analogous argument 
gives 

W^e(ClC2)A(«)(Cl) = A(^^)(Cl)W.(C2). (2.66) 

Alternatively one can note that the two equations (2.6) are the transpose of each 
other. However, in the non-symmetric case of section 6, this will relate the action of 
Ane and A^w on the vertex operators. 




Figure 2.7: Mapping from Ramond to Neveu-Schwarz sector using W^{C). 
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Now consider the products Ve{C)We{C') and We{C)Ve{C') which stand on the left 
hand side of figure 2.8. Using the star-triangle relation in the form shown in figure 2.3, 
we may move couplings H or M repeatedly from right to left, thereby interchanging 
the Ising couplings, and with them the rapidity variables. Because of the fixed spin 
boundary conditions this will eventually cease to have any further effect beyond mul- 
tiplication by the simple factor e^. Thus we have "commutation relations", although 
it is the rapidities rather than the operators which are interchanged. This is required 
by the fact that they intertwine the Neveu-Schwarz and Ramond sectors. 

Consider first the product Fe(C,^') = Ve{C)We{C)i we have 

i?e^(^/^')F,(C, CO = J2^''^^^''^-'fAC', C), (2.7a) 

e' 

where R comes from the first application of the star-triangle relation. 

We(C') We'(0 




We'(Q We'(0 
Figure 2.8: Commutation of vertex operators. 

Typically we shall require combinations of the form F+{(,(') ± F-{(,('), for which 
we have 

i?e^(^/^')(F+(C,0 + i^-(C,C')) = coshH{C/C){F+{C,C)+F-{CX)) 
i2e^(^/^')(F+(C,0 - i^-(C,0) = sinhiy(c/0(i^+(C',C) - i^-(C',C)) 

When we are dealing with expectation values rather than operators, we shall take 
quotients of such factors, as a result of which R will play no further part. From the 
star-triangle relation one sees that 

tanhif(C) = tanhL(C), 

the formula for which was given in (2.2). 

The case of the product -F'ee'(C)C') = ^e(C)K'(C') is similar, and we find 

e^(^/^V,,,(C,C') = e^(^/^')-'F,,KC',C)- (2.76) 
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For this formula tanhM((') = tanhi^(C). We shall, however, only use this relation in 
the case that e = e', for which we have the simple result that the rapidities may be 
interchanged without introducing any factors. 

3. Bulk magnetization 

As a preliminary to the calculation of surface and corner magnetizations in later 
sections, we will develop the method, and the underlying assumptions, by deriving 
the Onsager result for the bulk magnetization of the Ising model. 




Figure 3.1: Lattices with a defect line. 



3.1 Basic properties. Consider a lattice with a defect line as shown in figure 3.1. The 
pictures represent the unnormalised expectation values 

F,(Ci,C2) = trace(A(«)(a:)V,(C2)A(^^)(x)Vr,(Ci)), 

Fi(Ci,C2) = trace(A(^5)(a;)M^,(C2)A(«)(x)V;(Ci)), 

where we have combined pairs of CTMs, using the group property (2.5), to eliminate 
the bulk rapidity. The relation between F and F' is a simple reflection about the 
horizontal axis. This interchanges the couplings K and L, which is equivalent to 

c — ^ C 5 

Fe(Cl,C2) = i^e'(^/Cl,^/C2). 

This can be obtained directly from the definitions since trace(X') = trace(X) for any 
matrix X. 

It is easy to see that the dependence on and ^2 is only through the ratio Ci/C2- 
This is because these expectation values are matrix elements of a spin operator be- 
tween asymptotically degenerate maximal eigenvectors of the horizontal row transfer 
matrices above and below the spin. But, because of the star-triangle relation, all 
eigenvectors of the row transfer matrices except the line of dislocation are indepen- 
dent of rapidity while the latter can only depend on the ratio C1/C2 [12]. For the 
homogeneous case the magnetization is 

° \F+{c,o+F-{c,oJ V n (c, + (c, c) ; ' 
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Therefore we define the function G(Ci/C2) as 



'"^'^''^'^ - U+(O.C.) + F-«i.«J • 

similarly for G'(Ci/C2)) and 

Mo = G{l) = G'{l). (3.2) 
Finally, we shall require an inversion relation, namely 

GiOGi-0 = 1, (3.3) 

from which the normalisation of G{Q may be determined. It follows from 

Fe(-Ci,C2)=peF,(Ci,C2), 



where p = i^^-nL g^j^^ g^j^g ^]^g number of bonds of type K and L in the 

vertex operator VFe(Ci), regarded now as finite but large. To derive it, note that the 
transformation Qi — > —C,i induces K — > K+i7r/2, L — > L—i7r/2 in the operator We(Ci) 
of figure 3.1, and the Boltzmann weights transform as exp(±i^) — > ±(z) exp(±i^), 
exp(±L) ±(— i) exp(±L). F^{±(i,(2) are sums over configurations of the system. 
For the ground state term the relation is obvious; after replacing (^i by — Ci, P accounts 
for the phase factors (±i) and the sign changes if e = — 1 since one pair of spins are 
not aligned in that case. But all other terms in the sum differ from the ground state 
only by having an even number of further spin pairs of W^((i) which are not aligned; 
so every term in satisfies the same relation. From that, the result (3.3) follows. 




Figure 3.2: Steps in deriving equation (3.4). 



3.2 Functional equation. Deriving the functional (^-difference) equation involves quite 
a few steps which we describe in detail. Start with the function F^{(i,(2), and perform 
the following steps. 
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(i) Rotate the vertex operator V'e(Ci) 180° anticlockwise, using the boost property- 
together with the cychc property of the trace, to get the left hand picture of 
figure 3.2. As a result, (i — > x^i. 

(ii) Interchange the rapidities using (2.76), to give equality with the right hand pic- 
ture. 

(iii) Rotate W^{x(i) 180° anticlockwise to get the lattice which defines F^{x'^(i, (2)- 
As a result we find 

i^i(Cl,C2)-i^e(^'Cl,C2). (3.4) 

For the second relation we proceed as in figure 3.3: 

(i) Starting with -Fe(Ci, C2), rotate the vertex operator W^{(i) 180° anticlockwise to 
get the left hand picture of figure 3.3. Again we have (i — > x(i. 

(ii) Interchange the rapidities, this time using (2.7a), to get to the the right hand 
picture. This requires that wc use the coupling H{C2/xCi). 

(iii) Rotate V^{xC,i) 180° anticlockwise to recover the lattice which defines F^{x'^(^i, (^2)- 
Putting all of this together, we get the functional equations 

i?Fe(Ci,C2) = $]e^(^^/-^^)"'F;,(a;2Ci,C2). (3.5) 

e' 

Finally, by substituting either of (3.4, 3.5) into the other, we get uncoupled functional 
equations for either of the two functions. 




Figure 3.3: Steps in deriving equation (3.5). 



Changing to the functions G{$,), G'{$,) removes the dependence on the factor R 
which comes from the star-triangle relation, and gives functional equations in a single 
variable. We also set x = q^^'^ to get 



G{0^tanhH{l/qy'OG{q'0, 
G'(0 = tanhiy(l/g3/2^)G'(g2^). 



Ising model magnetizations 



11 



3.3 Solution of functional equation. We shall only discuss the solution for G(^). Using 
the formula (2.2), we obtain the basic functional equation, 



G(0 i-C/q; q'){-qy^; q^m q'W/C; q^) ' 
For the purpose of solving functional equations, note that the functions 

fiO = 9(0 = (3.6a) 

satisfy 

Here we employ the double infinite products 

oo 

iz;pi,P2)= n (i-^pW)- 

ni,n2=0 

In the present case this gives 

(qC; q\ q')hqVC; q\ q'){-q\, 9^ q')iqVC; 9^ q') 



{-q^; q'^, 9^)(g^/C; q^, q^){q^^; q^, q^){-q^/^; q^, q^) ' 



Notice that a possible normalising factor depending on k is fixed (to unity) by the 
inversion relation G{^)G{—^) = 1 of equation (3.3). Substituting now into (3.2), we 
arrive at the formula 

Mo = Gil) = 
(g^;g^,g^)(-g;g^,g^)' 

where half the terms already cancelled. Noting further that the remaining ratios fit 
into the patterns of (3.6), we arrive at Onsager's original result 

Mo = j^^^ = (1 - k')'/\ (3.7) 

So, although the form (1 — fc^)^/^ is simple, the infinite product carries a message 
about the imdcrlying structure of the problem. It is interesting to note that Yang 
derived precisely this infinite product in his original paper in which he proved the 
Onsager formula [2]. 

3.4 Uniqueness of solution. The fundamental assumption of the method is that the 
function G(^) is meromorphic in ^ and that its analytic properties are determined 
by the analytic properties of the Ising couplings as parametrized in (2.1). That is, it 
is a doubly periodic function of u. Even so, the functional equations by themselves 
determine their solution only to within a "pseudo-constant"; that is, a solution of 
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G{^) = G{q^^). Such a function is doubly periodic in the variable u so it is fully- 
determined by the structure of its poles and zeros. 

Our basic hypothesis is that the analytic structure of G{^), which is actually the 
quotient of two functions F±, is completely determined by the factor tanhif which 
occurs in the equation for that quotient. In this case the solutions will only consist 
of the infinite product of those factors, which also will automatically impose the 
normalisation up to a constant depending only on k. This constant was determined 
from the inversion relation (3.3). The calculations, and the arguments employed, are 
exactly parallel to those of [15] for the spontaneous polarization of the eight- vertex 
model. This approach also has some precedent, although not so explicitly stated, in 
Baxter's book [12]. For example, the working of pp231-236, to obtain the free energy 
of the eight- vertex model, is exactly the solution of a functional equation (10.8.24) 
of the type we have just solved, and the solution is the infinite product of just those 
terms which determine the relation. 

In the case just cited [12], arguments are presented that the solution actually has 
the appropriate meromorphic properties. Actually, for the Ising model it is possible 
to construct an explicit analytic expression for G{^) as an integral [25], using the 
fermionic description of the CTMs [21,22], and check that it does have the correct 
structure. 

For the results which we present in the remainder of this paper we know of 
no such integral expressions, but there are many points of contact with previous 
results, some analytic, many numerical, and in all cases we have carefully verified 
that there is exact agreement. It is also the case that all of our results can be viewed 
as special cases of a single function. But the result for a (10) surface has long been 
known [3], moreover one of our results for a (10) corner has recently been obtained 
independently by solving integral equations [8] . Each of these checks is for the values 
of the meromorphic function obtained here along a smooth curve in the complex 
plane, which is very strong evidence for the validity of our approach. 

4. Surface magnetizations 

In this section we generalize the method to a lattice with a (11) surface. This 
will allow us to calculate the magnetizations in the first and second surface row. 

4.1 Boundary reflection. The essence of the bulk calculation is a rotation of vertex 
operators through a total of 360°. For a system with a boundary this is not possible, 
but instead one can use a procedure in which the vertex operators are rotated back 
and forth, being "refiected" from the boundary after a change of the rapidity. This 
mechanism forms the basis of all further calculations. 

We are interested in simple free boundaries and we now show how this case can 
be handled in an elementary way. Consider the system shown in figure 4.1. Free 
boundaries require that we sum over the edge spins. The effect of this is manifested 
through terms of the form 




cr 
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where cr is a boundary spin, and cri, (72 the two second row spins joined directly 
to it by K and L. Now the crucial point is that the right hand side is completely 
symmetric with respect to K, L and also ui, a2- The interchange of all pairs K, L in 
a vertex operator which is adjacent to a free edge is effected by ( — > x/(; this gives 
the following simple construction: 

a vertex operator adjacent to a free edge can have its rapidity "reflected" off 
the boundary via the transformations Ve{C) ^ Veix/C), We{C) We{x/(). 
This avoids the complication of having alternating rapidities and "double row transfer 
matrices" which are normally associated with the use of boundary reflections and K 
matrices [26]. 




Figure 4.1: System with a free boundary. 

From this reflection property follow immediately two symmetries of the expecta- 
tion value shown in figure 4.1: 

FeiCi, C2) = F,{x/Ci, C2) = F,(Ci, X/C2). (4.1) 

Finally, we note that the inversion relation (3.3) still holds, since the ground state 
configuration, even with free boundaries, still has all spins except e fixed to +, and 
the argument of section 3 then goes through. 

4.2 Functional equations. We now derive functional equations for the system shown 
in figure 4.1. In order that we can re-use these equations for a corner, we temporarily 
replace x by an arbitrary constant p in the CTM A^^^\x/Q, dealing instead with 
A^^^\p/Q. The boundary reflection property, however, still involves the variable x. 

We shall also consider directly the function G(Ci,C2), as for the bulk magneti- 
zation. For the first functional equation, we require the following steps (see figure 
4.2): 

(i) Reflect the rapidity Ci off the right boundary; — > x/Ci- 

(ii) Rotate the right vertex operator 180° anti-clockwise so that the two are adjacent; 

(iii) Interchange the rapidities by multiplying G with the flrst factor tanhi?(Ci^2/pa;). 

(iv) Reflect the rapidity px/(i off the left boundary; px/d Ci/P- 

(v) Interchange a second time multiplying with the second factor tanhiy(^i/p^2)- 

(vi) Rotate back 180° to get the original system with (i Ci/p^- 
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(4.2a) 



(iv) 



We(C2) 
Ve(Ci/p) 



(ii) 




Ap/Q 




(iii) 

Ve(px/Ci) 

(V) 



We(Ci/p) 



Aq 



Ap/Q 




— o 
" Hi ^ 



Ap/Q 




H2 Hi 



Figure 4.2: Steps in deriving the functional equation for a free surface. 



A similar sequence of moves involving the other vertex operator gives 
G'(Ci,C2) = tanhif(x/KiC2)tanhif(Ci/K2)G'(Ci,p'C2)- 



(4.26) 



4.3 Solution of the functional equations. We have quite a number of relations which 
must be solved in a consistent way. The first point to notice is that there are two 
factors which determine the solution of the functional equations; one depends on the 
ratio C1/C2 and the other on the product CiC2- Therefore we use the Ansatz 

G'(Cl,C2) = *(Cl/C2)*(ClC2). (4.3) 

Imposing the symmetries (4.1) on this Ansatz gives 

$(1/0 = $(0, *(^/0 = *(^0, *(0 = H^/x), (4.4) 

so we need only find a single function ^{^). Equations (4.2) give two functional 
equations for each of $(0 and "^(C), however it is easy to check that they are all 
equivalent to the single equation 



$(e) = tanhiy(e/p)$(e/p'). 

Using (3.6) and requiring that $(0$(— = I5 obtain 



(4.5) 



^ ^ (pg^/y^; p^ Q^) (-PQ^/^^; p^ Q^) (-PQ^^V^; P^ q^) jpq^^'^C; g^) g^ 

(_pgl/2/^. p2^ g2) (p^3/2^. p2^ (pg3/2/^. p2^ g^) (-pg^/^C; P'^ , Q'^) 
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Recall that, in order to use these results for an edge, we must set 

p = X = g^/^. 

4.4 First row magnetization. The result for G(Ci, C2) gives directly the magnetization 
in the second row, which will be discussed below. But it also contains the result for 
the edge, because upon rotating Ve and We into a vertical position, as shown in figure 
4.3, the spin e appears in the first row. For a homogeneous lattice, the rapidities then 
have to equal x/^. Their original values before the rotation thus must be = 1) 
(2 — X, as indicated also in the figure. As a result, the edge magnetization is 

Ml = $(l)$(a;). 




Ve(x) We(l) 



Figure 4.3: Magnetization in first row. 

For both functions the double infinite products can be simplified into single infinite 
products as in section 3, giving 



{-q; g2)2(_gi/2.g)- 



According to (3.7) the first factor equals (1 — A;^)^/^. The second one differs from it 
by the substitution q q^^"^, which corresponds to a Landen transformation k — > 
2/c^/V(l + k) of the modulus [27]. This leads to 

Ml = (1 - k^/^, (4.7) 

which corresponds exactly with the result in [28] , which was obtained for the isotropic 
case. We see, however, that Mi is independent of the anisotropy, as it must be, since it 
is the matrix element of a spin operator between the boundary state and the maximal 
eigenvector of the row transfer matrix, both of which are independent of This 
result also shows directly the surface exponent Ps = 1/2. 

4.5 Second row magnetization. This is obtained by setting Ci = ^2 = C ^(Ci) C2)) 
so that 

M2 = $(i)$(cVa^), 
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which does depend on the anisotropy. The dependence can, however, be simplified 
very much. By reducing the double products to single products using (3.6), then 
further to theta functions by (2.3), one first shows that 

■/./^l/2 S^(^(^ + /V2)) 



^C/x) = -i{k') 



cn(i(w + /72)) 



Using addition theorems for elliptic functions and substituting the Ising couplings 
from (2.1) finally leads to 

Ms = (l-/c)^/^coth(K + L). (4.8) 

This implies a simple relation between M2 and Mi which can also be derived directly 
[29]. The expression (4.8) is symmetric in K and L, and M2 exceeds Mi for all finite 
positive K, L, as it should. Only in the anisotropic limit, where one diagonal line of 
couplings becomes dominant, one has M2 — > Mi. This can already be seen from the 
initial expression for M2 by setting C — 1 and using $(a;) = $(l/a;). 
In the isotropic case (^^ = x) one has 

Mi'° = (1 - A;2)i/2 = (4.9) 

which was already noted in [28]. An analogous relation between bulk and surface 
quantities was found in [20] for the eight-vertex model. In that case the edge polar- 
ization (in the isotropic case) is the square of the bulk polarization. One can show 
that this result implies (4.9) at the decoupling point. 

The dependence of M2 on the anisotropy is physically not obvious, but it is 
exactly what one expects from an eigenvector calculation, for which we would calculate 
the matrix element between the maximal eigenvector of the row transfer matrix and 
the vector obtained by acting on the boundary state with the row transfer matrix. 
However, such a calculation would be technically quite intricate compared with the 
present method. 

5. Corner magnetization 

In this section we consider 90° corners formed by (11) surfaces and calculate the 
magnetization for three different spin positions. 




Figure 5.1: Truncated corner with free boundaries. 
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5.1 Truncated corner. We first consider the situation, shown in figure 5.1, namely a 
corner to which two vertex operators are attached, so that the tip is missing. This is 
the analogue of the edge geometry of figure 4.1, and all considerations of the previous 
section can be taken over immediately. The only difference is that now the vertex 
operators are rotated only through one CTM, so that in deriving the formula for 
G^(Ci)C2) one must set 

P = C 

The boundary reflection property, however, still involves the variable x. Therefore 
the equations (4.2,4.3,4.4,4.6) are unchanged and the magnetization of the centre 
spin e has the form 

M2(C) = $(i)$(q^/VC'). (5.1) 

We postpone a discussion of the anisotropy dependence and first look at the isotropic 
case, C = x^/^ = q^l^. Then 

^iso ^ (Q3/4.gl/2^g2)4(_^7/4.^1/2^^2)4 

2 (_g3/4.gl/2^g2)4(^7/4.gl/2^g2)4 

_ (^3/4.^2^4(^5/4.^2)4 

- (^3/4.^2)4(^5/4.^2)4 

The products can be expressed as a ratio of Jacobian elliptic functions with the help 
of (2.3), 

(g3/4.g2)(g5/4.g2) ^ _ /^sn(3ir/4) 
(_^3/4. ^2) (_^5/4. ^2) cn{m'/A) ■ 

Using Landen transformations to obtain the special values of sn and cn finally leads 
to the result 





Figure 5.2: Neveu-Schwarz and Ramond corners. 



5.2 Ncveu-Schwarz corner. This corner, shown on the left of figure 5.2, is obtained by 
adding a single spin with the coupling K to the truncated corner. The magnetization 
M2 is not changed by this extra coupling while the magnetization Mi = (e) is given 

by 

Mi(C)=tanhi^(C)-M2(C). (5.3) 
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Again, this simplifies very much for the isotropic case C = x^^'^ = q^/'^. One 
can then either reduce the products introduced by tanhK(C) as before, or simply 
substitute tanhK{Q = {y/k + ^/l + k)~^ to obtain 

5.3 Ramond corner. This geometry is shown on the right of figure 5.2, where it is 
also indicated that one can obtain it from figure 5.1 by rotating VFe(Ci) through the 
quadrant. Therefore the magnetization Ms of the spin e follows by putting (i — 1, 
C2 = C in G'(Ci,C2) so that 

M3(c) = mnc/x), 

The isotropic case is readily extracted since the double infinite products appear in a 
form which enables their reduction to single products by (3.6), leading to 

M|«° = (l-/c). (5.5) 

This formula was already given in [28]. By contrast, the other analytical results are 
new. In particular those for Ml^° and M2^° have been verified against the numerical 
calculations carried out previously in [28,30]. 

5.4 Critical behaviour. From the results (5.2, 5.4, 5.5) one reads off directly the corner 
exponent /?c = 1 of the isotropic system. However, for the corner considered here, 
Pc is not a constant but depends on the anisotropy (as it depends in general on the 
opening angle of the corner). This dependence can be calculated in a straightforward 
way from the exact result (5.1). While there is no general simplification of the double 
infinite products in that case, it is not difficult to find the asymptotic form near the 
critical point (A; — > 1) and hence evaluate Pc- 

The function $(C) consists of factors which have the form 

0(C V) - '^^'^^'^^^'^"^'^^'^^^ 



For k ^ 1 one has [27] I ^ 00, I' ^ n/2, so that q as well as and therefore the 
quantities p, ^, rj, all approach one. Therefore we set 

q = e-% P = e-"% ^ = e-2^% r; = e-^^'S 

with e ~ — 7r^/ln(l — k) ^ for k ^ 1, whereupon 



, , ( tanhim + na + x)e 

in0(C,r;) = 2^ In <^ 

m,n=0 ^ 



tanh(m + na + y)e 



For the asymptotic form we replace the double sum with a double integral and apply 
the approximation 

tanh(s + A) _ ^ 2A 
tanh(s) sinh2s' 
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Then 



n0(4,?7)~ — JJ^ sinh2(s + t) ~ 8e V a J 



More generally we see that 



n 



8ae 

i=i i=i 



^ln0(^i,?7j) ^^i^i-Vi)- 
Applying this to we have 

«(C)«exp(-^). (5.6) 

which is dependent on p via a, but not explicitly dependent on the argument C- 
Setting finally p = C fi^^s a = u/tt. With the value of e given above it follows that 
the corner exponent is 

Since one is near k = 1, the parameter u can be expressed from (2.1) as 

where Lc is the critical coupling across the corner. This is precisely the result one 
obtains by rescaling the anisotropic system onto an isotropic one with effective opening 
angle 9 = 2u, together with the conformal prediction fic = /'20 for the isotropic case 
[6,7]. The present direct calculation is the first analytical proof of this argument. 
Because $(C) in (5-6) does not depend explicitly on the same result for j3c follows 
from M3, as it should. 

Finally, a treatment of corners with other opening angles would just lead to other 
values of p. For example, the straight surface corresponds to p = x = q^^'^^ so that 
q; = 1/2 and one recovers the surface exponent f3s = 1/2. Thus (5.7) is the general 
result and u indeed plays the role of an effective angle. 

5.5 Magnetization curves. In order to give an impression how the corner magnetiza- 
tion varies with temperature for different values of the anisotropy, we have calculated 
Ml numerically from equation (5.3). The results are shown in figure 5.3 for five values 
of the ratio L/K of the couplings (rather than (). One can see that for large L/K the 
magnetization curve rises only slowly as the temperature is decreased. This mirrors 
the almost one-dimensional character of a system with strong couplings across the 
corner. Conversely, Mi increases rapidly for small L/K. The variation of the critical 
behaviour as described by Pc is just a particular aspect of this general pattern. 
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Figure 5.3: Magnetization Mi as a function of temperature 
for different values of anisotropy L/K. 



6. Surfaces in the principal direction 

The two previous sections have treated surfaces and corners formed by edges 
taken in the diagonal (11) direction. In this section we shall show how to treat edges 
in the principal (10), (01) directions using a layered system. 




P4 P3 P4 P3 Pi P2 Pi P2 

Figure 6.1: Layered system with a free boundary. 

6.1 Layered system. Consider the system shown in figure 6.1. It is essentially the 
same as figure 4.1, except that we introduce vertical rapidities Pi which alternate 
between columns. The meaning of the CTMs and vertex operators which appear in 
the figure are exactly as before, except that the couplings are calculated using the 
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ratio of horizontal to vertical rapidities, viz 

K = K{C/p), Ce{C,Ci,C2}, 

L^L{(/p), P ^ {Pl,P2,P3,P4}■ 
We are particularly interested in some special choices of the p's which produce 
(10) and (01) surfaces. These are shown in figures 6.2 and 6.3. Consider for example 
figure 6.2(a). To get it from the general layered system we must choose pi and p^ 
so that the couplings which belong to their columns take the limits K{(/p) — > oo, 
L{C,/p) — ^ 0. For this we must set pi = Pa = Q/x. This effectively removes half 
of the spins in the lattice, since each pair which becomes identified by our choice is 
then equivalent to a single spin. Then we choose P2 = = 1 since we do not wish 
to alter the couplings in the other two columns. Finally we must set Ci = C2 = C 
to get the homogeneous lattice with the (10) surface shown in figure 6.2(a). These 
specializations will be made after the general solution has been determined. 




Figure 6.2: Surfaces with (10) and (01) free boundaries. 

Similar considerations give figure 6.2(b). So for figure 6.2 we need. 

P2=Pz = 1, Pi=Pi = C/x, (6.1a) 

Pi=Pa = 1, P2=Pz = C> (6.1fe) 




Figure 6.3: Corners with (10) and (01) free boundaries. 

Figure 6.3 shows two types of corner, 90° and 270°. For these the appropriate values 
are 



P2=Pa^ 1, Pi = C/x, P3 = C, (6.1c) 
Pi=P3 = 1, P2 = C, P4 = C/x, (Q.ld) 
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We write A{(;ps,p4), B{(;pi,p2), We{Ci;pi,P2), V;(C2;P3,P4) for the operators 
shown in figure 6.1. Just as in section 2, these operators are related by crossing 
symmetry. Remembering that the adjoint reverses the direction of transfer, and that 
the crossing symmetry must now be written as Q/p ^ xp/ we see that 

b{c-p,p') = a'{x/cMpMp% 
v,{C;p,p') = w',{x/c-MpMp'). 

We need the boost properties and commutation relations for these operators. A short 
derivation is given in the appendix; here we state just the properties used in this 
section. Since we only have to rotate the operator W^{C,i;pi,p2) through 180°, we 
require the composition of two boosts. It is 

= w-eKi;P3,P4)^(^^nC;P3,P4)s(^^)(C;pi,P2). 

Next, we need the commutation relations (2.7) for the vertex operators which 
have now aquired extra rapidities pi. The important point is that the 180° rotation 
changes the vertical rapidities pi, p2 to pa, p4, so they play no further role in the 
commutation because they will match for the two operators which must be multi- 
plied; this means that equations (2.7) actually remain valid with the same choice of 
couplings, which depend only on the horizontal rapidities. 

The last ingredient is the "boundary reflection". Recall that this is simply the 
interchange of adjacent pairs K{(/p), L[C,/p') which are joined to a common spin on 
the free boundary, since the sum is symmetric in both the variables and the other 
pair of spins. This means that we want to make the replacements C/P ^ ^p'/Cj 
C/p' ~^ xp/C- At first sight this brings the difficulty that the vertical rapidities 
have changed place. But we observe that the replacement may also be written as 
C/p ~^ {^PP' / /P-i c/p' ~^ {p^PP' IO)Ip' 1 which is equivalent to a change in only the 
horizontal rapidity by C — > {xpp' /Q. Thus we have the simple rule: 

a vertex operator adjacent to a free edge can have its rapidity "reflected" 
off the boundary via the transformations V^(C;Pi)P2) ^ V^{xpiP2/C',Pi,P2), 

W,{C,P3,P4) ^ Weixp3P4/C-:P3,P'l)- 

6.2 Functional equations. We now derive functional equations for the system shown 
in figure 6.1. We shall re- use the figures and steps which were used in section 4, except 
that we retain the original meaning of x = g^/^ in the CTMs. For the first functional 
equation, we require the following steps (see figure 4.2): 

(i) Reflect the rapidity Ci off the right boundary: — > xpip2/Ci- 

(ii) Rotate the right vertex operator 180° anti-clockwise so that the two are adjacent: 

XP1P2/C1 x'^PiP2/Cl- 

(iii) Interchange the rapidities by multiplying G with the factor taiih.H{(^i(^2/x^pip2)- 

(iv) Reflect the rapidity x^piP2/Ci off the left boundary: .x^piP2/Ci ~^ P3P4Ci/^PiP2- 

(v) Interchange a second time multiplying with the factor ta,nh H (p^p^d / XP1P2C2) ■ 
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(vi) Rotate back 180° to get the original system with — > C,i/r'^^ where 

2 2P1P2 

r = X . 

P3P4 

From these steps we obtain 

G(Ci,C2) = tanhiy(CiC2/a;Vp2)tanhiy(a;Ci/r'C2)G'(Ci/r-',C2). (6.3a) 
A similar sequence of moves involving the other vertex operator gives 

^(Ci, C2) = tanhiy(p3P4/CiC2) tanhi?(a;Ci/r2C2)G'(Ci, r\2). (6.36) 

These equations, and equations (6.5) below, generalize those of section 4, to which 
they reduce if all Pi = 1. 

6.3 Solution of the equations. We follow the procedure of section 4 by using the 
Ansatz 

G'(Cl,C2)=*(Cl/C2)*(ClC2). (6.4) 

The boundary reflection property imposes symmetries which are a generalisation of 
(4.1), namely G(Ci,C2) = G{xpip2/Ci,C'2) = G{Cii XP3P4:/C2), which manifest them- 
selves as 

*(1/P3P40 = *(PlP20> ^{xpiP2/0 = ^(XP3P40, *(0 = ^{xpm/O- (6-5) 

As before we need only find a single function $(^), since equations (6.3) are all 
equivalent to 

$(0 = tanhiy(a;e/r2)$(e/r2). 

The solution may be written down by reference to equations (4.5,4.6), whereupon 
one sees that it is only required to make the substitutions p — > r, ^ — > x^/r = q^^^^/r 
in order to obtain 

(-rVe;r2,92)(^2^.^2^^2)(^^2/^.^2^^2)(_^^.^2^^2)- ^ • ) 

The actual solution for any particular surface or corner is obtained by substituting 
this result, together with the special values of the parameters pi, into (6.4, 6.5). 

6.4 90° corner. The relevant parameters for this corner are given in (6.1c) from which 
r'^ = X = q^/'^. For the magnetisation, we have 

Me(C) = $(l)$(l/C), r = q'/\ (6.7) 

In the isotropic case, ( ~ 9^^^, we have already evaluated both products; $(1) in 
section 4.4 and ^{q~^/'^) in section 5.1. Reusing those results, we have 



$(1) = Vi^, Hq-'/^) = ( ^— ^ ) (v^ + VT+k) , 



1/2 



1 + Vk 
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giving 

Mi'" = {l-Vk){Vk + VTTk). (6.8) 

We can show that the general formula (6.7) agrees with the known result of 
Kaiser and Peschel [8, 10], 

Mc = l-(cothi^-l)(cothL-l)/2. (6.9) 

Using the simple identity (cothK — 1)^ = 2(cosh2i^ — sinh2i^)/(cosh2i^ — 1), it may- 
be written in terms of half-angles as 

MJC) = (1 - k)— — , \ , \ -r-^. 7 — r^- (6.10a) 

'^^^^ ^ ^dn(m/2)(dn(m/2) +iA;sn(m/2)cn(m/2)) ^ ^ 

To see that this is the same as (6.7), first note that when we set = q^^'^ in $(1/C), 
the double products simplify to elliptic functions as follows: 

^ (9^/\; gV^ g^)(-gVC; g^)(-g^/\; g^)(g/C; g^) 

^'^^ i-q'/H; q'/', q'){qyC; q^l\ q^W^K: 9^/^ q^){-qK\ q'/^ q") 
iq'/\;q')iqC;q'){q/C;q'){q'/VC;q') 
i-q'/H; q')i-qC; q')i-q/C, q^)i-q'/yC, q') 
_ , sn{i{u/2 + r/2)) sn(z(w/2 + /')) 
cn(z(w/2 + r/2)) cn(z(w/2 + /')) " 

Using addition theorems to eliminate the quarter period /' gives the form 

M iC) = {l- k) ^^(^^/^^ dn(zt./2) - z(l + k) sn{iu/2) 

"^^^ ^ ^dn(m/2)(cn(m/2) -isn(m/2)dn(m/2)) ' ' 

After a fair amount of further calculation, one can show that the two formulae 
(6.10a, 6) agree. 

6.5 Surface. Using (6.1a) we have now = a;^ = q, and for the magnetisation, 

M,(C) = $(1)$(1/C), r = gi/2. (6.11a) 



Substituting r'^ = q in the formula for $ and reducing the double products to single 
products gives 

$(1/^) = iqCjq'MC-^ f) ^ k''/^nd{iu/2). (6.116) 



It is easy to reconcile equations (6.11a, b) with the formula of McCoy and Wu [3] 
for the magnetization at a (10) surface, namely 

1/2 

, cosh 2K — cosh 2L \ /„ 
^'=1 cosh2i^-l I ■ ("l^) 
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where L is the dual couphng; cosh2L = coth2L. One simply substitutes from (2.1), 
and uses half angle formulae for the Jacobian elliptic functions, to obtain 

= (1 - A;2)i/2nd(m/2). 

Finally, we note that choosing = (2 = C different from the bulk rapidity C 
leads to the system shown in figure 6.4. By summing out the spins in the first row 
one then obtains a (10) surface with modified couplings in the boundary, a problem 
treated in detail by Au-Yang [31]. 



Figure 6.4: System with a modified (10) boundary. 



6.6 270° corner. Using (6. Id) we have now = = q^^'^, and for the magnetisation. 

Me = $(l)$(g^/VC), r = q^/\ 

However, there does not appear to be any simple reduction to Jacobian elliptic func- 
tions in this case. This is evident from the fact that the two bases in the doubly 
infinite products, g"^/^ and q"^, are not related by Landen transformation in the simple 
way they were for the 90° and 180° cases. However, the asymptotic analysis of the 
section 5.4 is easily applied since the function $(^) of this section is essentially the 
same as there, and one sees that the critical exponent is 

Pc = 1/3, 

independent of the anisotropy, as it must be. 
7. Conclusion 

We have presented an approach which permits the calculation of order parameters 
in Ising models with various kinds of simple boundaries. The method uses CTMs 
which already on geometrical grounds are the natural tool for the shapes studied. 
But while a previous attempt to use them was only partly successful [30] , the vertex 
operator technique leads to a very elegant and compact solution. Moreover it has 
a marked geometrical character which makes it easy to use. We have deliberately 
presented it in some detail, because so far the method is not widely known. By 
contrast, we have not addressed questions of the infinite lattice limit which is always 
involved implicitly when CTMs are introduced. 
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The results obtained from the functional equations all have the form of infinite 
products which reflect the spectrum of the CTMs as well as the matrix elements with 
the boundary state vectors. In general the formulae are quite complicated. But they 
are not difficult to derive, and we have seen that amazing simplifications may be 
achieved which can lead to very simple final expressions for the order parameters. It 
can in fact be quite difficult to find these reductions, and in a way, when found, the 
simple forms obscure the origin and structure of the results. 

One can think of many applications to the Ising model beyond those studied 
in this paper. Here we have not studied two-point correlation functions, or indeed 
more general correlation functions, although they are solutions of similar (but more 
complex) functional equations. Other surface configurations could also be studied: 
for example, by adding more vertex operators at the boundaries one can in principle 
treat spins in deeper layers. Again, the layered system in section 6 has only been used 
for simple special cases related to a particular problem, and certainly not investigated 
in any generality. In this context one should remember that CTMs were themselves 
developed by considering inhomogeneous systems, as were vertex operators and their 
functional equations. Now, in this paper, inhomogeneties have been used to rotate 
the surface from the (11) to the (10) and (01) directions. 

The method of q'-difference equations is not confined to the Ising model, having 
been used in papers on the eight vertex model [15, 32], the Zn+i x Z^+i generalized 
Heisenberg model [33], and for the XYZ chain with a boundary [20]. However the 
method is presented here in the language of spins for the first time, and this simplifies 
matters considerably compared with using ABF models, or taking the eight-vertex 
model at the decoupling point. As such it applies naturally to the chiral Potts model, 
although there are formidable difficulties in that case due to the fact that the rapidities 
lie on a high-genus Riemann surface. 
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Appendix 

Group property. We need to extend the group property (2.5) to include alternating 
vertical rapidities. The derivation is the same for both Neveu-Schwarz and Ramond 
sectors, and for the NW and NE corners. Here we choose A^-^^\(;p,p') and show 
that 

where A^^^\Q — A'^^'^''((^; 1, 1). A simple derivation follows from the observation 
that A^^^\^;p,p') may be regarded as the infinite product of vertical semi-infinite 
row transfer matrices [12], 



A(^^)(Ci;p,/) = • ■■W,{xp'/Ci)V,{xp/Ci)W,{xp'/Ci)V,{xp/Ci). 
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Multiplying on the right by A^^^\(2) and invoking (2.6) repeatedly (to commute the 
V and W through A) we arrive at the desired result. One can write it also with a 
scalar normalisation factor, but we assume that such factors are absorbed into the 
definition of the CTMs, as for the original group property (2.5). 

Boost property. We sketch the derivation of the boost property (6.2). As in section 
2, it is trivial to see that 

^e(C;Pi,P2)^(^^)(C;Pi,P2) = A(^\C;Pi,P2mi\ 

;pi,P2) = A'^^^\C;pi,p2)Weil). 

We have also the original boost property (2.6) for Ve{C) and W^{Q, which may be 
used, together with (Al), to "lift" {A2) to 

K(CiC2;Pi,P2)A(^^^(Ci;Pi,P2) - a(^)(Ci;pi,P2)K(C2), 

W^e(ClC2;Pl,P2)A(«)(Cl ■,Pl,P2) = A^^^\Ci;pi,p2)We{C2)- 

Alternatively, on can view this as an expression of Z-invariance [17] applied to the 
rotation of a rapidity line through 90°; in this process the alternating vertical rapidities 
present in the horizontal row are absorbed into the CTM. 

Taking the transpose of {A2) gives the analogous property for B^^^\(^;p.,p'): 
S(^^)(Ci;pi,P2)Vr,(CiC2;Pi,P2) = W^e(a:C2)S(«)(Ci;pi,P2), 

B^''HCi;Pi,P2)V,{CiC2;Pi,P2) = v;«2)S^^^nCi;Pi,P2). 

Notice the factor x which comes from the crossing symmetry by which we obtain the 
;p,p') corner from the A^^^\(^;p,p'). The required result (6.2) now follows 
from the composition of (^2) with {A3). 
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